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Abstract. We study the three-pion decays of the lightest isoscalar vector mesons, lj and <j), in a dispersive 
framework that allows for a consistent description of final-state interactions between all three pions. Our 
results are solely dependent on the phenomenological input for the pion-pion P-wave scattering phase shift. 
We predict the Dalitz plot distributions for both decays and compare our findings to recent measurements 
of the (^ — >■ Stt Dalitz plot by the KLOE and CMD-2 collaborations. Dalitz plot parameters for future 
precision measurements of a; — > Stt are predicted. We also calculate the mr P-wave inelasticity contribution 
from iUTT intermediate states. 



PACS. 11.55.Fv Dispersion relations - 
interactions 



1 Introduction 



. The treatment of meson decays into three pions using dis- 
' persion relations is a classic subject. It was developed first 
in the 1960s in the context oi K -^ Sir decays [1], and 
already in the 1970s, it was applied specifically to the 
I decay w — > Stt [2]. One of its main virtues is the fact 
. that final-state interactions among the three pions are 
' fully taken into account, in contrast to perturbative, field- 
, theory-based approaches; the constraints coming from an- 
' alyticity and unitarity are respected exactly to all orders. 
This becomes the more important, the higher the mass 
. of the decaying particle, hence the higher the possible en- 
ergies of the two-pion subsystems within the Dalitz plot: 
while decays like X — >• Stt or 77 — > Stt have successfully 
been analyzed in perturbative settings like chiral pertur- 
bation theory [3, 4] or even non-relativistic effective theo- 
ries [5-7] , it is obvious that these are doomed to fail for the 
decays of w or into three pions, where the influence of 
the p resonance is already strongly felt (a;) or even shows 
up in the form of resonance bands inside the Dalitz plot 

In recent years, there has been a flourish of attempts 
to also treat the physics of vector mesons in effective- field- 
theory approaches [8-15], motivated not least because of 
their prominent nature in hadronic physics involving vir- 
tual photons, their relatively low masses among the var- 
ious meson resonances, and their strong coupling to the 
light pseudoscalars. However, when it comes to the three- 
pion decays of w and 4>, descriptions in terms of vector- 
meson dominance and improved tree-level models (mostly 
incorporating a finite width of the intermediate p res- 
onances) [13, 16-18] (see also Ref. [19] for earlier refer- 
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ences), which do not fully respect analyticity and unitarity 
constraints, still seem to be state of the art.^ 

In this respect, we consider this an ideal time to 
take up dispersive studies of these decays once more. 
High-precision phase-shift analyses of pion-pion scatter- 
ing are now available [20-22] and can be employed as 
input for decay studies. On the experimental side, high- 
statistics Dalitz plot investigations have either been per- 
formed {(j) — >■ 37r [23,24]), or are planned or ongoing 
(w — >■ Stt [25,26]). Finally, these systems are ideal test 
cases to study the achievable precision of Dalitz plot de- 
scriptions with theoretically rigorous methods: they are 
simple in terms of partial waves, as entirely P-wave domi- 
nated, and kinematically in the transition regime between 
the light (kaon, eta) decays and those of heavy {D, B) 
mesons that might be treated with similar methods in the 
future [27]. 

The outline of this article is as follows. We introduce 
the necessary basics on kinematics and partial-wave de- 
composition in Sect. 2, and describe the formalism for a 
dispersion-theoretical description of the decays w, cj) —> Sir 
in Sect. 3. First numerical results of our solutions are dis- 
cussed in Sect. 4, before we compare in detail to the exper- 
imental Dalitz plot studies for (j> —>■ in and predict Dalitz 
plot parameters for w — >■ 37r in Sect. 5. As a further appli- 
cation, we calculate the contribution of wtt intermediate 
states to the inelasticity in the pion-pion P partial wave 
in Sect. 6. We conclude in Sect. 7. Several aspects that 



^ It needs to be pointed out, though, that Lagrangian-based 
approaches often have the advantage of relating various differ- 
ent processes to each other; such symmetry constraints most 
of the time have nothing to do with analyticity and unitarity, 
and hence can at best be imposed a posteriori in dispersive 
studies. 
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lie somewhat outside the main line of argument are rele- 
gated to the Appendices, where we discuss higher partial 
waves (Appendix A), the analytic properties of the disper- 
sive integral (Appendix B), as well as possible contribu- 
tions from heavier resonances (Appendix C), and derive 
a generalization of the Omnes representation taking into 
account inelasticities (Appendix D). 

Preliminary results of this study have already been 
presented in Ref. [28]. 



2 Kinematics and partial-wave decomposition 

In this work we consider the threc-pion decay of the light- 
est isoscalar vector mesons, 

y(pv)^7r+(p+)7r-(p_)7r"(po) , V ^ uo / c^ , (1) 

where the particle momenta are conventionally related to 



the Mandelstam variables by s 
V+f, u = {pv - P-f, with 



{Pv -Po) , t = (pi. 



3so^s + t + u = M^ + 3M^ . (2) 

Here and in the following we restrict ourselves to the 
isospin limit, M^^o — M^± = M^r, unless explicitly stated 
differently. 

Since F — )■ Stt is of odd intrinsic parity, the amplitude 
can be decomposed as 

M (s, t, u) = ie^i,ai3n^P+P°:Po^{s, t, u) , (3) 

where n'^ is the polarization vector of the decaying vector 
particle, and T(s,t,u) is a scalar function. The absolute 
value squared of the amplitude reads 



\M{s,t,u)\^' 



1 r 



f2N2 



(4) 



—K^s)smHs\T{s,t,u)\'' 
16 



where 9s is the center-of-mass scattering angle in the s- 
channel, cos 0s = {t — u)/k{s), and 

K(s) = a,(s)Ai/2(M2,M2,s) , (5) 

with the Kallen function \{x,y,z) ^ x^ + y^ + z^ — 2{xy + 
yz + xz) and o-^(.s) = ^1 - AW^/s. 

Due to Bose symmetry only partial waves of odd angu- 
lar momentum can contribute to the amplitude. Accord- 
ingly, the partial- wave decomposition of the scalar part of 
the amplitude reads [29] 



T{s,t,u) = ^ fe{s)P^{zs) 



(6) 



£odd 



where Zs — cos6's, Pg{zs) is the differentiated Legendre 
polynomial, and one can project onto the partial-wave am- 
plitude fi{s) by 



1 f^ 

fds) = 7^ / dzs[Pe-i{zs) - Pe+i{zs)] J^{s,t,u) 



(7) 



For the dominant £ ~ 1 partial wave we have 



/i(s) = -y dz.,{l-z^,)^is,t,u). (8) 

In our analysis we neglect discontinuities of F- and higher 
partial waves (see the corresponding discussion in Ap- 
pendix A), so that T{s,t,u) can be decomposed in terms 
of functions of a single variable as 



J"(s, t, u) = T{s) + T{t) + F{u) , 



(9) 



where J^(s) only has a right-hand cut. The dispersion re- 
lations derived in the following will be given in terms of 
these single- variable functions J^{s). Equation (9) repre- 
sents a symmetrized partial-wave expansion, its general- 
ized form allowing for F-wave discontinuities is shown in 
Appendix A, Eq. (48).^ 

We will also make use of the partial- wave expansion of 
the TTTT -^ TTTT Scattering amplitude, which is convention- 
ally given as 



r^(s,t,u) = 327rV(2£+l)4(s)P,(z) , (10) 



•E 

i=0 



where i^ is the partial-wave amplitude of isospin / and 
angular momentum i and can be expressed in terms of 
the phase shift Sj according to 



tiis) = 



2i(T7r(s) 



(11) 



3 Dispersive analysis of three-body decays 

In the following, we will present our framework to account 
for three-body effects in a non-perturbative fashion. It re- 
lies on the formalism laid out in the context of dispersive 
studies of ?7 — > Stt using the same method [32, 33]. 

Setting up dispersion relations for three-body decays 
is a non-trivial task due to the very nature of the pro- 
cess in question: three-body final states necessitate the 
incorporation of three-particle cuts, a substantial analytic 
ingredient that is left out in isobaric models, like vector- 
meson dominance models [13,16]. The main idea behind 
the approach presented here is to set up a set of dis- 
persion relations for the corresponding scattering process, 
i.e. Vtt ^ TTTT, with My < 3M^ and s > {My + M^f , 
and analytically continue the resulting expressions to the 
physical realm of the decay process. My > ZMt, and 
AMl<s<{Mv-M^f. 

We wish to remark that we essentially only consider 
elastic unitarity in the following (see Sect. 4.4 for an 
attempt to partially account for inelasticity effects in 



^ Similar amplitude decompositions, often quoted as being 
due to a "reconstruction theorem" , have been demonstrated in 
the perturbative context of chiral perturbation theory to two 
loops for pion-pion scattering [30,31], and subsequently been 
extended to 77 —> Sir decays [32] . 
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Fig. 1. Contribution of KK intermediate states to — >■ 3n. 
The general rescattering KK — >■ 7r'''7r~7r'^, denoted by the thick 
dot in the left diagram, is entirely dominated by the reso- 
nance for {pk + Pk)^ = -^Ji S66 the right diagram. 



the pion-pion P partial wave). One might suspect that 
the dominant decay (p — > KK with subsequent inelastic 
rescattering KK — )• Stt, see Fig. 1, may have a significant 
impact on the Dalitz plot distribution for (^ — >■ Stt. This is 
not the case for the following reason: the KK intermediate 
state occurs at a fixed total energy of {pk +PkY = -^^Ji 
where the rescattering KK — > Stt will be entirely dom- 
inated by the (very narrow) resonance. The diagram 
therefore factorizes, and the whole effect of the KK in- 
termediate state can be absorbed into a complex wave- 
function renormalization constant for the field, which 
amounts to an unobservable overall (constant) phase of 
the amplitude. The Dalitz plot distribution of —> Stt re- 
mains therefore entirely unaffected. [A similar argument 
was already put forward to disregard inelastic rescattering 
contributions in 77' — > "q-Ki: in Ref. [34].] 

Note finally that in the sense of a dispersion-theory 
analysis, there is no pi: two-body intermediate state_ con- 
tributing to the (j) decay in a similar fashion as KK de- 
picted in Fig. 1: pTT is no distinct state from Stt, the effect 
conventionally encoded that way shows up in our analysis 
as the resonant two-body tttt P-wave interaction. From the 
point of view of dispersion theory, there is no meaningful 
way to differentiate between pir and Stt final states in </> 
decays. 



3.1 Unitarity relation 

In the following, we relate the discontinuity of the Vtt — > 
TTTT s-channel partial-wave amplitude to the amplitude it- 
self and the tttt P-wave phase shift. Similar considerations 
for t- and u-channel proceed completely analogously. In 
the elastic approximation with only pion-pion intermedi- 
ate states one has for the discontinuity of the diagram in 
Fig. 2 [35] 



disc7W(s, Zs) = - 



{2n) 



.M{s,z'JT'*{s,z':) 



x{2n)d{l'~M^){27r)S{{q-ir~M^), 

(12) 

where q — p+ + p-, z^ — cos ds, z[. = cos d'^ , and z'J = 
cos 9", where 9s denotes the center-of-mass scattering an- 
gle between the initial- and final-state momenta, 9'^ be- 
tween initial and intermediate state, and 0" between in- 
termediate and final state. We remark that in the case of 
a three-body decay (as opposed to elastic two-body scat- 
tering), the discontinuity across the cut is in general not 




Fig. 2. Diagram of the Vtt — >■ tttt s-channel scattering ampli- 
tude. The double line denotes the vector meson, all single lines 
refer to pions. The gray circle corresponds to the Vn — >■ tttt 
amplitude M{s, z'J), the white one to the tttt — > tttt amplitude 
T^{s,z's) (of isospin I = 1). To calculate the amplitude one 
cuts through the intermediate pion propagators with momenta 
I and q — I. 



purely imaginary, so that the usual relation discA^(s) = 
2iImA^(s) does not hold. To avoid confusion, we will ex- 
clusively refer to the discontinuity of the amplitude. 

Upon insertion of Eq. (3) and elimination of the d- 
functions, we can rewrite Eq. (12) according to 



diacT{s,Zs)^^^ I dz'^ I del)', \_y 



-1 Jo 



xT{s,zi)T'*is,z':), (13) 



where z'J — Zgz'^ + y^l — z1\/\ — z'^ cos0(,. Finally, we 
can use the partial-wave expansion for the amplitudes T 
and T^ given in Eqs. (6) and (10), and project out the 
i — \ partial wave to arrive at the Watson-like unitarity 
relation, 

disc/i(s) = 2ih{s)9{s - 4A//2)sin5(s)e"'''(') , (14) 

where (5(s) = 5\{,s) is the tttt P-wave phase shift. Noting 
that disc /i(s) — disc J^(s), the previous expression can be 
recast into 

disc J'(s) = 2i(j^(s)+ J^(s))6'(s-4Af^) sin ,5(s)e-*''(") , 

(15) 
where T(^s) is referred to as the inhomogeneity of the in- 
tegral equation and is given by angular averages over T 
according to 

^(.) = 3((l-z2)-F)(s), 
{z"f)(s)^-j dzz'7(5(3so-s + z«(5))) . (16) 

J^(s) contains the left-hand-cut contribution to the partial 
wave /i(s). 

The previous expressions are well-defined in the scat- 
tering regime. However, special care has to be taken in the 
decay region due to the non-trivial behavior of the func- 
tion k{s). It imposes a complex analytic structure on the 
angular integration, which is discussed in more detail in 
Appendix B. Aside from that Eqs. (15) and (16) form an 
integral equation, the solution of which we discuss in the 
following. 
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3.2 Homogeneous equation, pion vector form factor 

Before discussing the solution of Eq. (15) in its entirety, 
let us consider an obvious simplification, the homogeneous 
problem with J^(s) = 0. The unitarity relation then closely 
resembles that of the pion vector form factor F^ (s), and 
Eq. (14) is a manifestation of Watson's final-state theo- 
rem [36] : the form factor shares the phase of the (elastic) 
scattering amplitude. The solution to this problem is given 
in terms of the Omncs function [37] , 



FY(s)^p(s)n{s), 



f2{s) 



s 
cxp< - 



ds 



, Sis') 



TT y4jv^2 S'{S' - S) 



where P{s) is a polynomial and the Onines function is 
normalized to 12(0) — 1. From Eq. (17) we observe that 
the Omnes function is completely determined by the phase 
shift. In particular, the asymptotic behavior of the Omnes 
function is constrained by the one of the phase shift: for 



S{s — > oo) — >■ kiT, one finds f2{s) 



We shall assume 



k — 1 for the asymptotic behavior of the P wave, which 
guarantees the high-energy fall-off for the form factor ex 
1/s as suggested by perturbative QCD (up to logarithmic 
corrections) exactly if P{s) is a constant, and hence due 
to gauge invariance, P{s) — 1.^ 

Pertaining to this issue the behavior of the tttt phase 
shift is not known to arbitrarily high energies. In our nu- 
merical calculations we use a phenomenological parame- 
terization of the phase shift up to a certain energy s — Ajj, 
beyond which it is set to a constant. At the accuracy at 
which we are working this constant does not have to be ir 
exactly. This issue is studied in Sect. 4.4. 

One possibility to suppress the high-energy behavior 
of the phase shift in the Omnes function is to resort to 
a twice-subtracted dispersion integral, as opposed to the 
once-subtracted version used in Eq. (17) (cf. e.g. Ref. [41]). 
Due to the identification F^ (s) = ^(s), the additional 
subtraction constant can be related to the charge radius 
of the pion. 



n{s) 



cxpi g(r2)];'s+'* 



Sis') 



ds' 

I" J4M^ S (■?' - s) 



(18) 



Comparing to Eq. (17), we can express the charge radius 
in terms of a sum rule. 



2\V 



{rX 



6 



ds 



,5is') 



I" J4M2 



(19) 



To take advantage of the suppression of high energies in 
the oversubtracted dispersion integral, one may make use 
of an independent phenomenological determination of the 
charge radius. This allows us to reparameterize our lack of 



^ Note that for a precision analysis of the pion vector form 
factor up to about 1 GeV one needs to account for p-u> mixing 
and the onset of inelasticities. This is not of interest in the 
context of our work and we refer to Refs. [38-40] for further 
details. 



knowledge of the behavior of the amplitude at large s as 
a polynomial, which should be a decent approximation at 
the energies we are working at. In Sect. 4.4 we discuss the 
numerical effects of using the two different versions of the 
Omnes function. We want to emphasize that, as a mat- 
ter of principle, using a radius different from the sum-rule 
value Eq. (19) is inconsistent: it leads to a wrong (expo- 
nential) high-energy behavior. In practice, however, and 
for the (small) deviations in the charge radius we con- 
sider, the effects of this error are not visible within the 
physical region of the decays under consideration. 



(17) 3.3 Integral equation and solution strategy 



To find a unique solution to the unitarity relation Eq. (15), 
one derives an integral equation for J"(s)/.f?(s) instead of 
Tis) (the issue of ambiguous solutions to the unitarity 
relation is sketched in Refs. [32,42]). One finds 



J"(s) sa\Sis)Pis) 
disc -ttA = , „; ,, ' ti[s 



fii& 



\nis)\ 



AMi 



(20) 



the solution of which is readily obtained using Cauchy's 
integral formula. 



J"(s) ^[2is){a + 



ds' sinSis')Tis') 
TrJ,^,.~\Qis')\is'-s)\ ' 



(21) 



where a is a subtraction constant. The order of the sub- 
traction polynomial is limited by the asymptotic behav- 
ior of the integrand. The Froissart bound [43] constrains 
the behavior of the amplitude for large s, Al(s,t, m) < 
Cs log (s) for some constant C, and consequently -F(s) < 
C"s~^/^ log^(s) for s — > oo (and another constant C"). 
Along with the asymptotic behavior of the Omnes func- 
tion it is obvious that the integral remains finite. We ob- 
serve that since !F is linear in F we can further simplify 
Eq. (21) from the point of view of the numerical imple- 
mentation: 



J"(s) = aFais) , J"(s) = aFais) , 

s f°" ds' sin Sis')Fais') 



Fais) = Qis)il + 



^J,Mi s' \Qis')\is' - s) 



(22) 



This observation is extremely useful, since our results will 
be pure predictions aside from an overall normalization 
constant that can be fixed after the iteration process. For 
lack of a better theoretical method to fix the normaliza- 
tion, we fit to the experimentally determined partial decay 
width unless explicitly stated otherwise. 

The integral equation (22) can now be solved by an 
iterative numerical procedure: we start from an arbitrary 
input function Fis) and calculate the inhomogcneity Fis) 
from Eq. (16), which in turn is used as an input for the 
calculation of an updated Fis) by means of Eq. (22). The 
process is repeated until the solution converges to a fixed 
point with sufficient accuracy. The final result is indepen- 
dent of the specific choice of the starting function. We will 
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use J-'{s) ~ i7(s) as our starting point: as the tttt P wave is 
dominated by the p resonance at low energies, this closely 
corresponds to an isobaric description of the decay, and 
the modification of J^{s) in the iteration procedure, the 
difference between starting and fixed point, allows us to 
quantify crossed-channel effects generated by the iteration 
in a plausible way. 



3.4 Oversubtraction 

The integral in the solution of the dispersion relation for 
y — > Stt decays, Eq. (22), is guaranteed to converge, given 
our assumptions on the high-energy behavior of ampli- 
tudes and phases. This solution has the maximal degree 
of predictability, as it only depends on one single real pa- 
rameter, the subtraction constant a that merely represents 
the overall normalization of the amplitude (the phase of 
which is of course unobservable) ; the complete Dalitz plot 
distribution is then a prediction. 

We will discuss various sources of the theoretical er- 
rors in this representation, most of which are in one way 
or the other associated with the high-energy behavior, in 
some detail in Sect. 4.4. It is obvious, though, that while 
our high-energy constraints on the amplitudes, including 
e.g. the assumed smooth continuation of the scattering 
phase, present a plausible and internally consistent proce- 
dure to interpolate between the very well-constrained low- 
energy part and the asymptotic behavior as suggested by 
the Froissart bound, we certainly neglect various details in 
the description of an intermediate-energy range, in partic- 
ular inelastic contributions. The hope (which, eventually, 
has to be checked phenomenologically) is that inside the 
dispersive integrals, this intermediate-energy range does 
not influence the low-energy description of the decay am- 
plitudes too much. However, similar to what we explained 
for the Omnes function, in order to suppress the influ- 
ence of inelastic contributions even further, we can alter- 
natively subtract the dispersive solution once more than 
strictly necessary, at the expense of introducing another 
subtraction constant: 



T(s) = n(s)\a + h's + 



ds' smS(s')J'{s') 

This is equivalent to Eq. (21) if b' obeys the sum rule 
,, 1 p ds' sin S{s')J'{s') 



"^ J4M2 S ■ 



Ms')\ 



(23) 
(24) 



Because of the special analytic structure of J^(s) discussed 
in Appendix B.2, which is due to three-particle cuts in the 
decay amplitude, the subtraction constant b' is complex. 
If one allows it to take values different from the sum rule 
in order to give the dispersive representation more free- 
dom in a fit to experimental data, it therefore represents 
two new real parameters, modulus and phase of b' . Again, 
as in the discussion for the Omnes function, such a b' dif- 
ferent from its sum-rule value is in principle inconsistent 



and leads to a high-energy behavior violating the Froissart 
bound. We will, however, again find that these violations 
do not manifest themselves in practice in the energy range 
considered here. 

The linearity of Eq. (23) in the subtraction constants 
a and b' also massively simplifies the numerical solution 
strategy in this case.'* The full solution can be constructed 
as the linear combination 

J^{s)=a[j^:{s)+bMs)] , (25) 

where b ~ b' /a, in terms of the basis solutions found from 
"2 r-oo ds' sin S{s')f'^{s')' 



J-^(s) = l2(.)|l + ^ 



4M3 s'^ \[2{s')\{s' ~ s) 



Tb{s)^f2{s)\s + 



2 /-oo II ^, 



ds' sin S{s')Tb{s') 
n J, Mi ^Ms')\{s'-s) 



(26) 



■F'ais) and J^h(s) can therefore be calculated by the it- 
erative procedure explained above, before adjusting the 
subtraction constants in a fit to experimental data. 



4 Numerical results 

In this section we show the numerical results from solving 
Eq. (22). We start off by discussing the numerical input 
in Sect. 4.1 before showing the convergence behavior for 
both a; — > 37r and (/) — > 37r in the iteration procedure in 
Sect. 4.2. In Sect. 4.3 Dalitz plot distributions are pre- 
sented along with a study of crossed-channel rescattering 
effects, followed by a discussion on how these effects hold 
up against the expected errors of our analysis in Sect. 4.4. 
The comparison to experiment follows in Sect. 5. 



4.1 Numerical input 

The integral equation Eq. (22) is fully determined except 
for the TTTT input and the subtraction constant. This in- 
put is subject to uncertainties, the discussion of which we 
defer to Sect. 4.4. Instead, we only give a central set of 
parameters used in the following. It was already pointed 
out that we use the partial decay width to fix the subtrac- 
tion constant (that serves as the overall normalization of 
the amplitude), namely [44] 



r„ 



uj^-Stt 



7.56 MeV , r^^g^r = 0.65 MeV 



(27) 



It should be noted that we do not consider errors on the 
partial decay widths, since the uncertainties thus gen- 
erated are by far superseded by the error sources dis- 
cussed in Sect. 4.4. Furthermore, the masses involved are 
given by M„ = 782.65 MeV, M^ = 1019.46 MeV, and 
M^ = 139.57 MeV. We use the tttt P-wave phase shift 
based on an ongoing Roy-equation analysis [21] (first as- 
pects of which have recently been published [45] ) . There 



■^ We are grateful to Gilberto Colangelo, Stefan Lanz, and 
Heiri Leutwyler for pointing this out to us. 
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Stt decays 



An 

A 

# subtractions in Q(s) 

phase-shift param. 

inelasticities 



1.3 GeV 
2.0 GeV 



Ref. [21] 



Table 1. Input to the analysis, see text for explanations. 



are other parameterizations of the P-wave phase shift of 
comparable accuracy available, see for example Ref. [22]. 
We will discuss the influence of the difference between the 
two in our error discussion. 

There are other error sources that relate to the upper 
limit of the dispersion integrals in Eqs. (17) and (22). Our 
procedure in the Omnes integral is as follows: the phase 
shift derived from the Roy-equation analysis is strictly 
known up to the validity limit of the Roy equations of 
y/s ~ 1.15 GeV. We use the somewhat extended phe- 
nomenological parameterization of Ref. [21] up to Aq = 
1.3 GeV, set the phase to a constant beyond that point, 
and calculate the Omnes integral analytically. The upper 
limit of the dispersion integral in Eq. (22) has less physical 
significance, it is rather an indicator how well one sums up 
the remainder of the integral. In our analysis the integral 
is cut off at yl = 2 GeV. We emphasize that this does not 
mean that we know the physics of tttt interactions up to 
that point, but it certainly is a better approximation to 
the integrand than setting it to zero. 

We have also considered methods to include clastic res- 
onances between 1.3 and 2.0 GcV (p'(1450) and p"(1700)), 
as well as inelasticities e.g. from Air intermediate states. 
Furthermore we have estimated the possible contribution 
of a p3-dominated tttt F wave. The corrections stemming 
from these contributions are tiny and also deferred to the 
error discussion and the Appendices. Our input parame- 
ters are summarized in Table 1. 



4.2 Convergence behavior of the amplitude 

To illustrate the convergence behavior of the iteration pro- 
cedure, the amplitudes T{s) for w — > 37r and </) — > Stt are 
plotted after each iteration step in Figs. 3 and 4. Note that 
we show the normalized amplitudes, i.e. a = 1. Conver- 
gence of the cj — > 37r amplitude is reached fast, with F{s) 
all but indistinguishable from the final result after two it- 
erations. The same holds true for (j) — )■ 37r, although one 
more iteration step is required. However, in both cases the 
difference between the final result and the starting point is 
significant: since our starting point is the Omnes function 
that resums rescattering between two pious only, this hints 
at sizeable crossed-channel effects in the decay region. 

We observe that the p peak in w — )■ 37r is slightly en- 
hanced and shifted toward lower energies. Even though 
the peak is not part of the physical decay region and 
this shift therefore not observable directly, the somewhat 
steeper rise should leave its imprint on the Dalitz plot. 
The — > 37r amplitude exhibits a similar shift of the p 
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Fig. 3. Successive iteration steps of real (upper panel) and 
imaginary (lower panel) part of the amplitude Ta{s) for ct) — > 
Stt. The vertical dashed lines denote the physical region of the 
decay. 



peak toward smaller energies. In contrast to the lo de- 
cay, the peak here is attenuated, which should certainly 
have an impact on the Dalitz plot distribution, since it is 
part of the physical region. We reiterate that the modifi- 
cations we observe due to three-particle rescattering lead 
to non-negligible effects on the amplitude. In the following 
section we shall discuss how our observations translate to 
the actual Dalitz plot distributions. 



4.3 Dalitz plot distributions and crossed-channel 
rescattering 

For a three-body decay, the Dalitz plot distribution is of 
particular interest. The Dalitz plot is a two-dimensional 
scatter plot in terms of two kinematic variables; we use 
the dimensionless variables 



t — u 
VSRv 



y = 



sq - s 
Rv 



(28) 
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where Ry = |My(My - 3M^). 

In Fig. 5 we show the Dalitz plot distribution divided 
by the P-wave phase-space factor given in Eq. (4) and 
normalized to 1 in the center of the Dalitz plot. Figure 5 
is thus a pure prediction: it is free from any input aside 
from the tttt P-wave phase shift, which is well-established 
up to at least 1.15 GeV, thus covering the entire physical 
range for both processes. 

The (x> — > Stt Dalitz plot exhibits a relatively smooth 
distribution, which rises from the center to its outer bor- 
ders with a maximum increase of roughly 20% with respect 
to the center. The available phase space is not sufficient 
to contain the p resonance. This behavior unambiguously 
fixes the sign of the leading slope parameter in a possible 
Dalitz plot parameter representation to be positive, see 
Sect. 5.3 for a detailed discussion and numerical results. 

The — > Stt Dalitz plot in contrast shows significantly 
more structure, since the physical region encompasses the 
p resonance: the resonance bands of the p^ and p are 
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Fig. 5. Dalitz plots for a; — >■ Sir (upper panel) and — >■ Stt 
(lower panel), normalized by the P-wave phase space. 



clearly visible in Fig. 5. From its center, the Dalitz plot 
distribution rises towards these bands with a maximum 
enhancement of roughly 60% at the peak, and then steeply 
falls off, showing almost complete depletion towards the 
outer corners. 

In order to illustrate effects of crossed-channcl rescat- 
tering on the Dalitz plot distribution, we study the phase- 
space corrected Dalitz plot after the iteration procedure 
(|J^fuiip), divided by the same quantity before the itera- 
tions, corresponding to J^ = 0, i.e. the sum of pure Omnes 
solutions (|-7^jr=ol^)- ^^ devise two approaches of fixing 
the subtraction constants. First we assume that the sub- 
traction constant is given by some independent method, 
and we are interested in what bearings the crossed-channel 
effects have on both the overall shape of the Dalitz plot 
and the partial decay width. The quantity |.?"fuiiP/|-^;r=ol^ 
is then independent of the specific choice of the subtrac- 
tion constant. For the sake of the argument we fix the 
subtraction constant from the experimental decay width 
for J^ — and then run the iteration procedure. Our sec- 
ond approach is to readjust the subtraction constant in 
such a way as to reproduce the experimental decay width 
in both cases, with and without crossed-channel rescat- 
tering effects included. The focus then lies exclusively on 
changes to the profile of the Dalitz plot distribution. 

The results of both approaches are shown in Fig. 6. 
By keeping the subtraction constant fixed we observe that 
crossed-channel rescattering enhances the cj — > Stt partial 
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Fig. 6. |7-fuii|Vl-^^=oP for 



3n (upper panel) and 



Stt (lower panel). The two difTerent scales next to the plots 
correspond to two different methods of fixing the subtraction 
constants: a is fixed to reproduce the decay width before the 
iteration only (left scales), or it is fixed to reproduce the decay 
width before and after the iteration (right scales). 



width by roughly 20%. This qualitative behavior trans- 
lates to the Dalitz plot, where the enhancement of the 
distribution ranges between 10% on the borders and 24% 
in the center. The — > Stt partial width in contrast sees 
a decrease by likewise roughly 20%. The decrease across 
the Dalitz plot is stronger, with 60% on the border, than 
in the center where it amounts to roughly 10%. Fixing the 
subtraction constant before and after the iteration proce- 
dure, the effects due to crossed-channel rescattering are 
to a large degree absorbed in the partial width in both 
decays. Indeed, the remaining effect on the a; — )■ Stt Dalitz 
plot amounts to an 8% decrease that is alleviated towards 
the central region, where one observes a slight increase of 
roughly 3%. The same qualitative behavior is observed in 
(j) — )■ Stt, only quantitatively stronger: the 50% suppression 
on the border is counterbalanced by a 20% enhancement 
in the center. The p bands are left unscathed by the iter- 
ation procedure. 

Overall, we find that crossed-channel rescattering 
leaves a significant imprint on the Dalitz plot distribution. 
However, before checking how our approach to those three- 
particle effects holds up against experimental scrutiny, we 
study whether the size of the effects even withstands the 
uncertainties of our input parameters in the following sec- 
tion. 





cj — > Stt 


[%] 
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-> Stt [%] 


# subtractions 


5.2 






4.5 


S{s) parameterization 


2.4 






2.6 


Ao 


0.6 






0.8 


S{s) — > n 


0.4 






0.6 


inelasticities 


0.2 






0.4 



Table 2. Relative variation of Im J^(s) in the p peak as per- 
centage of our standard input parameters (Table 1) for the 
respective error source for a; — > Stt (left) and — > Stt (right). 
See text for further details. 



4.4 Error discussion 

As we stated in the previous sections our results are — 
aside from the subtraction constant — fully constrained by 
the TTTT P-wave phase shift. The parameterizations of the 
phase shift we use are very accurate in the low-energy 
regime; in this section we study the influence of uncer- 
tainties that are mainly due to the high-energy behavior 
and different methods to assess it. These error sources are 
listed in Table 2 together with the maximum variation 
they produce in ImJ^(s) on the p peak relative to our 
central parameter set (see Table 1). In the following, we 
discuss these error sources in more detail. 

As mentioned in Sect. 3.2, using a twice-subtracted 
Omnes function is a means to suppress the high-energy 
behavior of the phase-shift input. If we naively plug our 
central parameter set into Eq. (19) we obtain 



{r\ 



V 

sum/ TT 



0.415 fm^ 



(29) 



which lies somewhat below a next-to-next-to-leading order 
chiral perturbation theory (ChPT) analysis [46], 



(-; 



ChPT/TT 



0.452 ± 0.013 fm^ 



(30) 



and the current particle-data-group average [44]. Disper- 
sive analyses of e+e^ -^ tt+tt^ data [39,40] point towards 
a value much closer to Eq. (29), with central values of 
the order of (r^)^ ~ 0.43 fm^; however, we consider the 
variation between a once-subtracted and twice-subtracted 
Omnes function, using the phenomenological radius in the 
latter, a conservative estimate of the uncertainty. We find 
that the variation of the p peak amounts to roughly 5% 
and is thus the largest source of uncertainty. 

It is essential to try to assess the uncertainties stem- 
ming from our main input, the parameterization of the 
TTTT P-wave phase shift. Since these phase-shift solutions in 
themselves are very accurate, especially in the low-energy 
regime, we address these errors by varying between two 
different parameterizations, referred to as "Bern" [21] and 
"Madrid" [22] parameterizations in the following. These 
lead to a difference in Im J^(s) in the p peak by roughly 2%, 
which amounts to the second largest error in our analysis. 
We have also checked how higher resonances (/9'(1450) and 
p"{1700)) can modify the phase shift (beyond the range 
of the Roy analyses) and hence our results. This aspect is 
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Fig. 7. Real (upper panel) and imaginary (lower panel) part 
of the amplitude J-a{s) for uj — > Stt. The vertical dashed lines 
denote the physical region of the decay. The shaded line repre- 
sents the uncertainty generated by upper and lower boundary 
of our solution (see text). 



described in greater detail in Appendix C; the impact on 
the decay amplitudes in the physical region however turns 
out to be negligible. 

Related to the previous point is the exact value of 
y/s beyond which the phase is set to a constant, Aq. 
To check the dependence on this error, we vary Aq — 
1.3 ... 1.5 GeV. Note that the phase shift beyond 1.3 GeV 
is not well constrained at all, yet it still produces a 
very small uncertainty below 1% and will be neglected in 
the following. An equally negligible modification is found 
when leading the phase continuously to tt beyond Aq in- 
stead of setting it to a constant, which is in principle in 
better agreement with the 1/s behavior of the tttt vector 
form factor for large s. 

Inelasticities in the tttt P wave have been incorporated 
in a simplified fashion, following the method proposed in 
Ref . [47] , by bestowing the exponential of the phase shift 



with an inelasticity parameter ri{s) == ?7i(*)i 



t\{s) 



2i 



(31) 



This leads to modified dispersion integrals (see Ref. [47] 
and Appendix D), 



Tis)=aas)S{s)f2{s) 



xU + 



ds' [e'^(^')-7?(50e-»^(^')]j-(gO 



(32) 



where 



rj ^/^(s) above the cut, 
C(s) ^ { 77^/2(3) below the cut, (33) 

1 elsewhere, 



and 



, log?7(s') 



-W=cxp|- f ds-^ 



l&Ml 



(34) 



where ^ denotes the principal value integral, and we as- 
sume that inelasticities set in at the four-pion thresh- 
old. The inelasticity starts showing major deviations from 
unity only above roughly 1 GeV (see also Sect. 6). Con- 
sequently the effects of using such a simplified model for 
the inelasticity show little impact (below half a percent) 
on our final result. 

We add one final remark on the integral cutoff in the 
dispersion integral of the full amplitude, Eq. (22). This 
cutoff has been fixed to 2 GeV and not varied. The rea- 
son for this is that it is merely a "numerical" cutoff: the 
kinematic range of validity is fixed by the tttt phase-shift 
parameterization. The dispersion integral has no physical 
content beyond that point; we just need to ensure that the 
integral cutoff is large enough so as not to produce numer- 
ical artifacts. We found that this is the case at 2 GeV. As 
a side remark, we have checked that lowering the integral 
cutoff as far as to 1.3 GeV still affects the amplitudes in 
the physical region by less than the main sources of un- 
certainty discussed above. 

After identifying the charge radius in the twice-sub- 
tracted Omnes function and the phase-shift parameteri- 
zation as the main single error sources, we will now per- 
form a more detailed study of their combined error. It 
turns out that we can identify certain parameter sets as 
upper and lower boundaries of the amplitude for that pur- 
pose. The combination of the Bern parameterization along 
with a once-subtracted Omnes function will serve as the 
lower boundary, whereas the upper boundary is given by 
the Madrid phase together with a twice-subtracted Omnes 
function and the determination of the pion charge radius 
from Ref. [46] . The uncertainty bands generated between 
these boundaries are shown for lj -^ Z-k and — >■ 37r in 
Figs. 7 and 8. The maximum uncertainty of the ampli- 
tude in the physical region amounts to roughly 8% in both 
cases. 
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Fig. 8. Real (upper panel) and imaginary (lower panel) part 
of the amplitude J-a{s) for (p — >■ Stt. The vertical dashed lines 
denote the physical region of the decay. The shaded line repre- 
sents the uncertainty generated by upper and lower boundary 
of our solution (see text). 



To rule out that interference of the crossed channels 
amplifies the uncertainties, we also analyze how the errors 
develop across the Dalitz plot. For that purpose, we study 
the following quantity across the Dalitz plot: 



l-^ul 



m\^ 



iJ-^p + iJ-zp 



(35) 



where J^u/i denotes the normalized amplitude on the up- 
per/lower boundary. We find that for a; — > Stt the error 
in the Dalitz plot is even slightly decreased with respect 
to the amplitude, the maximal uncertainties amount to 
1.4% towards the corners of the Dalitz plot. For ^- Stt 
we observe the largest uncertainties in the outer edges of 
the Dalitz plot, where the Dalitz plot strength itself is 
strongly suppressed. So while the relative error there rises 
up to 8%, the absolute error is very small. 

Comparing to our results in Sect. 4.3 we clearly see 
that crossed-channcl contributions are sizable enough to 



outweigh the uncertainties. In the following section we will 
analyze whether crossed-channel effects are actually ob- 
servable in experiment. 



5 Comparison to experiment 

5.1 Fit to the cf) ^ 3tv Dalitz plot: single subtraction 

We now wish to compare our theoretical predictions to 
the experimental results for the — > Stt Dalitz plot mea- 
surements by the KLOE [23] and CMD-2 [24] collabo- 
rations. The former has significantly larger statistics (al- 
most 2 X 10^ events in the Dalitz plot) than the latter 
(close to 8 X lO** events); furthermore, the energy reso- 
lution in the KLOE measurement was significantly bet- 
ter than the bin size (in the range 1-2 MeV compared to 
a bin size of 8.75 MeV), such that smearing effects were 
found to be negligible, and we could fit our amplitudes to 
efficiency-corrected data directly (with purely statistical 
errors based on data and Monte Carlo statistics), while 
for the comparison with the CMD-2 data, they had to be 
convoluted with efficiency matrices by the collaboration 
before. While consistency with both data sets is clearly 
desirable, we will present the comparison to the KLOE 
data in some more detail in the following. 

Our first goal is to perform a fit to the Dalitz plot 
distribution with our most predictive theoretical repre- 
sentation, Eq. (22), employing a single subtraction, such 
that the normalization is the only free parameter of the fit. 
The shape of the Dalitz plot is thus purely a prediction, 
and we can compare the x^ of the fit with and without 
crossed-channel rescattering. There are two caveats to this 
procedure, which we need to discuss beforehand. 

First, our calculations are performed in the isospin 
limit of equal charged and neutral pion masses; we use 
the charged pion mass, not least for consistency reasons 
due to the fact that the tttt phase shifts are only avail- 
able in the isospin limit, with the charged pion mass used 
as the reference quantity. The effect of this approxima- 
tion on the amplitude is expected to be small (compare 
e.g. Ref. [7] for a detailed study in the context of r/ — > 37r), 
with the main difference due to different charged and neu- 
tral pion masses showing up in the available phase space: 
the true physical Dalitz plot is slightly larger than in our 
calculation. To account for this dominant isospin-breaking 
correction, we therefore multiply the amplitude (squared) 
with the physical phase-space factor: 



|Al0^3^(s,t,M)|" 



16 



Kl{s)sm'^0a\T{s,t,u)\'^ , (36) 



where in contrast to Eq. (4) we have cos^o — it^u)/Ko{s) 
and 



Ko(s) 



/ 4 M^ 
'l-^^^Ai/2(M2,M2o,.) 



(37) 



Numerically, we employ M^o — 134.98 MeV. Furthermore, 
in order to avoid distortions due to threshold effects in 
the decay amplitude, we omit data bins that cross the 
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Ti 


T^ 


Ti 


:f^ 




(full) 




i^-- 


= 0) 


X^/ndof 


1.50 


1.17 


1.71 


2.06 


a X GeV^ 


228.6 ±0.2 


216.6 ±0.2 


200.2 ± 0.2 


187.6 ±0.2 


a„ X 10^ 


7.9 ±0.5 


7.4 ±0.5 


10.9 ±1.0 


12.6 ±1.2 


4>u 


-0.30 ±0.11 - 


-0.10 ±0.11 


0.87 ±0.06 


0.95 ±0.06 



Table 3. Fit results to the KLOE data for (f> — > Stt, using the once-subtracted dispersive representation. Shown are the outcomes 
of the fits for upper and lower boundaries of the theoretical uncertainty band, with and without crossed-channel rescattering 
included. The normalization of a is arbitrary, such that only the relative changes in a between different fits are significant, not 
the absolute values. The uncertainties quoted refer to the errors of the fit. 




Fig. 9. Selected data from the KLOE measurement [23]. 
Shown is the efflciency-corrected number of counts in the re- 
spective bin, divided by the phase-space factor in Eq. (36) and 
normalized to 1 in the Dalitz plot center. 



boundary of the Dalitz plot, see Fig. 9, thereby ensuring 
that our amplitude is never evaluated below the isospin- 
symmetric threshold, that is, for (M^± + M^aY < t, u < 
4M2±. 

The second caveat, discussed in Ref. [23], concerns the 
fact that the (f) is produced in e+e^ collisions at DA<?NE, 
e+e^ —>(/)—)■ TT+TT^TT^, which allows for the background 
process e^e~ — > wtt" — )■ 7r"'"7r~7r'', the lott'^ invariant mass 
equalling the mass of the (j), with the w subsequently 
decaying into tt+tt". As the decay a; — > tt+tt^ violates 
isospin, with the branching fraction suppressed to the per- 
cent level [44] , the overall modification of the Dalitz plot is 
small. ^ However, it is entirely concentrated in the narrow 
band s = M^, and leaves a visible effect there. There are 
two possible strategies to deal with this issue: one could 
simply omit the corresponding horizontal band in a fit of 
the Dalitz plot; or, alternatively, add a resonance term of 



"" Note that this is an isospin-violating effect specific for the 
37r production in e+e^ collisions. If we interpret the decay 
uj — >■ TT^TT^ in terms of a p-uj mixing angle 6pui (see e.g. 
Ref. [48]), this effect is linear in dp^^ due to the fact that the 
photon has both isoscalar and isovector components. In con- 
trast, isospin breaking due to p-co mixing in rnr scattering is 
necessarily suppressed to second order in 6puj and hence irrel- 
evant. 



the form (see Ref. [23]) 



M2 



M2 - i^r^ - s 



(38) 



to Eq. (9) after the iteration. (Note that we have factored 
out the normalization constant a for reasons of comparison 
with the fit results in Ref. [23].) M^^ and I^j are fixed to the 
particle-data-group values [44] , and a^^ and (j)^^ are dimen- 
sionless fitting parameters. We follow the latter strategy, 
in particular since the resonance term also has a small im- 
pact on bins adjacent to the horizontal band at 83.7 MeV. 
Our standard x^ fit is performed with 1834 data points 
and three free parameters. We perform separate fits for the 
upper and lower boundaries of the theoretical uncertainty 
band as discussed in Sect. 4.4, both with and without 
crossed-channel rescattering effects included. The results 
are listed in Table 3. The data set fitted to is given in ar- 
bitrary normalization, such that the constant a does not 
correspond to the real subtraction constant a and is only 
shown for comparison of the changes between different fits. 
We notice that the fit quality considerably improves once 
crossed-channel rescattering is taken into account. This in- 
dicates that crossed-channel effects lead to modifications 
that are not only non-negligible, but even observable in the 
structure of the Dalitz plot. It is also observed, however, 
that the precision of the data is such that the variation 
of the x^/ndof within the theoretical uncertainty is quite 
non-negligible: the upper boundary yields a considerably 
better fit than the lower boundary for the full crossed- 
channel analysis, while the opposite holds for the pure 
Omnes solutions. The a^^ coupling is found to be in the 
1% range, the expected order of magnitude; its numerical 
value determined together with the full dispersive solu- 
tion is in slightly better agreement with what is found in 
Ref. [23]. The best x^/ndof of 1.17 of our fits is not quite 
as good as it is in Ref. [23], and it needs to be pointed 
out that due to the high number of degrees of freedom, 
the p-value characterizing the goodness of the fit is still 
rather low even in the best case, p(x^ = 1.17) — 3x 10~^. 
However, this fit quality is achieved with less degrees of 
freedom: apart from the wtt" background term (which has 
nothing to do with the genuine — > 37r Dalitz plot), this is 
a one-parameter fit, the shape of the Dalitz plot in our case 
is a pure prediction. In Ref. [23] there are two additional 
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degrees of freedom by fitting the (complex) "background" 
term. From the construction of our dispersive ampHtude, 
it is obvious that in this form, such an independent back- 
ground term is inconsistent with the requirements of an- 
alyticity and unitarity. Furthermore, in Appendix C we 
show that a simphfied approach to including higher res- 
onances does not allow for an improvement of the fit as 
opposed to claims made for the nature of the background 
term being due to p' effects. It seems that crossed-channel 
effects saturate the background term in the KLOE data 
to a large degree. 

Figure 10 shows slices of constant y through the Dalitz 
plot, where the number of events per bin divided by the 
bin efficiency is plotted against the bin number (note that 
this is not the same as Fig. 5 in Ref. [23], where slices 
of constant x are shown). Our fit results are displayed 
as error bands. The full solution (white band with dashed 
boundaries) gives a better description of the data than the 
sum of three Omncs functions (gray band), particularly in 
the central region of the Dalitz plot. 

These conclusions on the significance of rescattering 
effects are not immediately substantiated by the compar- 
ison to the CMD-2 data [24]. The fits performed by the 
collaboration seem to lead to the almost opposite result: 
a fit based on the sum of Omnes functions leads to a very 
good x^/ndof of about 1.0 for both variants of the Omncs 
function discussed; while the full amplitudes yield bad fit 
qualities of x^/ndof = 1.5 . . . 1.8. We note, however, that 
due to the significantly smaller number of degrees of free- 
dom (ndof = 197), the p- value for the CMD-2 fit of the 
full Tu solution, p(x^ = 1.50) = 7x 10~^, is even better 
than the best KLOE fit in Table 3. Still, this is a sur- 
prising result in different respects: primarily, as the phe- 
nomenological fits of Breit-Wigner p resonance terms plus 
a constant "background" amplitude yield perfectly com- 
patible results in both experiments; secondly, as Ref. [24] 
cites the significance for a non- vanishing background term 
at 3.3(7, it is somewhat unexpected that simply replacing 
the Breit-Wigner p by an Omnes function is already suf- 
ficient to yield a good description of the data. 



5.2 Fit to the ^ ^^ Stt Dalitz plot: two subtractions 

In order to understand the situation of the two — > 37r 
data sets better, we attempt a description with the (more 
flexible) twice-subtracted representation Eqs. (25), (26). 
Due to the equivalence with the once-subtracted form 
in case the additional parameter h fulfils the sum rule 
Eq. (24), the fits can only improve: the hope is to find an 
acceptable fit also to the CMD-2 data that is compatible 
with the fundamental principles underlying the dispersive 
representation. We refrain from employing the two dif- 
ferent variants of the Omnes function defined earlier; we 
expect the second subtraction in the dispersion integral to 
render the second subtraction inside the Omncs function 
redundant, and therefore only use the standard form (17). 
We perform the fit of this representation with the two dif- 
ferent TTTT P-wave phase parameterizations [21, 22] to check 
the consistency of the outcome. 





Born [21] 


Madrid [22] 


X^/ndof 


1.02 


1.03 


a X GeV* 


207.6 ± 1.4 


207.1 ±1.5 


(sum rule) 


228.6 


225.4 


\h\ X GeV"'' 


0.97 ±0.03 


0.94 ± 0.03 


(sum rule) 


0.72 


0.75 


argb 


0.52 ± 0.03 


0.42 ± 0.03 


(sum rule) 


0.73 


0.70 


a^ X 10^ 


7.3 ±0.6 


7.5 ±0.6 


<i>ui 


0.40 ±0.10 


0.33 ±0.10 



Table 4. Fit results to the KLOE data for — > 37r, using 
the twice-subtracted dispersive representation. Shown are the 
outcomes of the fits using the phase parameterization of the 
Bern [21] and Madrid [22] groups. The "sum rule" entries for 
modulus and phase of h refer to the evaluation of Eq. (24) with 
the once-subtracted dispersive representation; the "sum rule" 
entry for a just serves as a reminder that also the normalization 
of the amplitude changes significantly between the two fits. 
For a, the same remark holds as in Table 3. The uncertainties 
quoted refer to the errors of the fit. 



The results of the fits to the KLOE data are shown in 
Table 4. The resulting values for a and 6 are compared to 
the ones obtained from the fit of the once-subtracted rep- 
resentation, see the previous section, via the sum rule (24) . 
For both phases, excellent fits of x^/ndof close to 1.0 are 
obtained (|5(x^) — 0.22 . . .0.24). Closer comparison shows 
that Im6 in fact stays very close to its sum-rule value 
(in particular for the Bern phase parameterization) , while 
Re h is shifted (enlarged) more significantly, at the expense 
of a somewhat reduced normalization a. We compare the 
resulting function Ta^s) + bTb{s) to the once-subtracted 
T{s) as well as the Omnes function in Fig. 11. We see that, 
indeed, the difference to the once-subtracted J-"(s) is small 
compared to the latter's error band, let alone the differ- 
ence to the Omncs function. We also display the fit result 
in comparison to the KLOE data in Fig. 10 as the black 
band. Note that as we now optimize the additional param- 
eter through the fit routine and only vary between the two 
different phase parameterizations, the band is significantly 
narrower than the other two. While the deviation from the 
once-subtracted fit results seems to be really minor, the 
overall improvement in the x^/ndof is significant. 

Fits to the CMD-2 data have again been performed by 
the collaboration. This time, they result in a x^/ndof of 
0.96 and 0.94 (with associated p-values of 0.64 .. .0.71), 
using the Bern and Madrid phases, finding 



.13 J "•^"J-0.12/ 



argfe== {0.00±0.16, --0.18±0.18} , 



(39) 



respectively. While the modulus of b therefore agrees per- 
fectly with the fit to the KLOE data, the phase seems 
to prefer a real h, in contradistinction to the sum-rule 
prediction. Only the KLOE data therefore yield a high 
significance for a non-vanishing phase of the parameter b. 
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Fig. 10. Efficiency-corrected number of events per bin plotted against bin number for slices of constant y. Every peak corresponds 
to one slice/one specific value of y, progressing from negative to positive x, and slices arc ordered from negative to positive y. 
Our fit results are plotted for the sum of Omnes functions J-" = (gray band), for the full solution of the once-subtracted ansatz 
(white band with dashed boundaries), and for the full solution in the twice-subtracted variant (black band), together with the 
KLOEdata [23]. 
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Fig. 11. Real (upper panel) and imaginary (lower panel) part 
of the amplitude J-{s) (normalized with a = 1) for — >■ Stt, 
comparing the fit of the twice-subtracted dispersive representa- 
tion to the KLOE data to the once-subtracted J-{s) as well as 
the Omnes solution. All curves are shown with the uncertainty 
bands attached as discussed in the text. 



We have not investigated systematic uncertainties in the 
determination of b from data. 



5.3 uj ^^ ^TT Dalitz plot parameterization 

Since there is currently no precise data available on the 
Lo ^ in Dalitz plot, we will now discuss the issue whether 
it is feasible to use a parameterization of the Dalitz plot in 
terms of a polynomial in a precision determination, and 
if so, with how many terms necessary. This is common 
practice for decays with final-state particles at low ener- 
gies, such as iT — > Stt, -q -^ Stt, or 77' -> rjmr. Since the 
Dalitz plot for uj -^ 7r+7r~7r° is relatively smooth, we ex- 
pect that a similar description should be possible here, 
even though the phase space available is somewhat larger. 
Moreover, as opposed to e.g. the decays K^ — > tt^tt^tt^'^. 



T] -^ Stt", and 77' -^ ?77r"7r°, one does not face the issue of 
cusp effects [5,6,34,49] in the decay region that give rise 
to possibly large non-analytic structures. 

The w — )■ Stt decay amplitude is fully symmetric un- 
der exchange of s, t, and u, see Eq. (9), such that the 
Dalitz plot description is formally similar to the one in 
1] — > Stt". For a parameterization in terms of a polyno- 
mial, it is therefore convenient to rewrite the Dalitz plot 
variables in Eq. (28) in polar coordinates, namely 



y 



' z sm ( 



X = -y/zcos^ . (40) 

We can therefore parameterize the amplitude squared in 
terms of z and (j) according to 

|j:-p„i(^, ^)|2 ^ |A/'|2|l + 2az + 2/3z3/2 sin3</> + 2-fz^ 

+ 2dz^^^sm3<j) + 0{z'^)\ , (41) 

where A^ is the normalization and a, /?, 7, S are the Dalitz 
plot parameters, in strict analogy to 77 — t- 37r° [7]. Note in 
particular that we attempt to parameterize |J^|^ in poly- 
nomial form, not |A^|^: the trivial kinematic factor due 
to the P-wave characteristic of the vector-meson decay 
should always be retained exactly. In the following, we 
study whether such a parameterization makes sense when 
trying to quantitatively describe the a; — > Stt Dalitz plot, 
and give predictions for the parameters. 

We perform a fit of the polynomial Dalitz plot rep- 
resentation \Tpo\{z , (j>)\'^ to our theoretically determined 
amplitude \Tthiz,(t')\'^j minimizing the function 



^' = i//^^ 



Af-D = / dzd(l} 

'■D 



3z(3so ~ -RwA/zsin3(/))\ 
(M^ + 3M,)2 ) 

|J-poi(z,0)|2-|J-th(^,0)|' 



|AA|^ 



(42) 



where T) denotes the area of the Dalitz plot. The normal- 
ization of Eq. (42) is chosen such that yx^ denotes the 
average deviation of the polynomial parameterization in 
|7W(s, i, 'u)|2, relative to the Dalitz plot center. The prefac- 
tor (in round brackets) corresponds to the kinematic fac- 
tor relating \M{s,t,u)\'^ to \F{s,t,u)\'^ in Eq. (4), rewrit- 
ten in terms of z and 0, and normalized to 1 in the center 
of the Dalitz plot. We decide to include the kinematic 
phase-space factor in the minimization in order to give 
less weight to the outer parts of the Dalitz plot, which are 
expected to also contribute statistically less in an experi- 
mental determination. 

We start with the singly subtracted dispersive repre- 
sentation, where the Dalitz plot shape is a full-fledged 
prediction. Our results are summarized in Table 5, where 
again we vary the amplitude between upper and lower 
error boundary of the once-subtracted solution, and also 
compare to the pure Omnes solution without rescatter- 
ing effects. We observe that for the full amplitude, a fit 
with two parameters (the normalization is not counted as 
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\Af\ X GeV^ 



ax lO'^ 



13 X 10^ 



7 X 10^ 



<5x 10-^ 



X^ X m' 



max dev.xlO 



full 


1451. 


.1447 


84... 96 (102) 


— 


— 


— 


0.9. 


.1.1 


40 




1453. 


.1449 


74 ... 84 (90) 


24 ... 28 (30) 


— 


— 


0.052 . 


.0.078 


3 




1453. 


.1449 


73 ... 81 (86) 


24 ... 28 (29) 


3... 6 (8) 


— 


0.038 . 


.0.047 


3 




1453. 


.1449 


74 ... 83 (88) 


21 ... 24 (25) 


0...2(3) 


7. ..8(9) 


0.012. 


.0.011 


2 


f = 


1433. 


.1429 


137... 148 


— 


— 


— 


1.1. 


.1.3 


50 




1435. 


.1431 


125 .. . 135 


29 ... 33 


— 


— 


0.25. 


.0.29 


20 




1436. 


.1433 


113... 120 


26... 29 


24 ... 27 


— 


0.036 . 


. 0.045 


4 




1436. 


.1433 


114... 122 


23... 25 


20... 23 


7. ..8 


0.002 . 


.0.003 


0.3 



Table 5. Fit results for the Dalitz plot parameters with (above) and without (below) crossed- channel effects included. We show 



the Dalitz plot parameters along with the values for \/x^ ^s defined in the text, "max dev." is the maximum deviation at any 
point across the DaUtz plot between the full solution and the polynomial fit. The numbers in brackets for the full solution refer 
to the extension of the parameter ranges considering the twice-subtracted dispersion relation; for details, see main text. 



a fit quantity) already gives a very good description of the 
theoretical data: compared to the one-parameter fit, the 



X^ is improved by roughly a factor of 15, the maximum 
relative deviation across the Dahtz plot, i.e. the quantity 
\J^poiiz,4')\^/\J^thiz,4i)\'^ - 1 at any point, is improved by 
approximately a factor of 10 from 4% to 0.3%. When in- 
creasing the number of Dahtz plot parameters beyond 2, 
the changes in the parameters themselves as well as the 
quality of the fit are rather small, and quite probably be- 
yond the reach even of a precision experiment. Note that 
the maximum deviation between two and three parame- 
ters is approximately the same, however the value of yx^ 
is significantly reduced. The reason for this is that the 
contributions with stronger deviation are pushed toward 
the outer boundary of the Dalitz plot, owing to the phase- 
space factor in Eq. (42). 

The Omnes solution shows a slightly different behav- 
ior: a good fit quality is reached only with three param- 
eters, and the third parameter 7 is found to be signifi- 
cantly larger. Deviations across the Dalitz plot with only 
two terms are still substantial at 2% and probably within 
reach of an experimental determination. Also the param- 
eters a and P still see changes in the 10% range. The 
comparison of our predictions of the Dalitz plot parame- 
ters between the full and the Omnes solution hints at a 
significant infiuence of crossed-channel rescattering. Note 
that, beyond the precise values of the Dalitz plot parame- 
ters, the sign of the leading parameter a is unanrbiguously 
fixed in a; — > Stt as discussed before (see Sect. 4.3). 

Given our comparison to the experimental (j) —^ 3tt 
Dalitz plots, we may still wonder how reliable these pre- 
dictions of w — > 37r Dalitz plot parameters are. It is ob- 
vious that as soon as we switch to the twice-subtracted 
dispersive representation, we cannot strictly predict all of 
these any more: at least a would have to be an input quan- 
tity. In order to estimate the potential effects, we resort 
to the following procedure: we assume the deviations of 
the second subtraction constant b from the sum-rule re- 
sult to be moderate; as these should be due to imperfectly 
understood high-energy behavior of the amplitudes, we 



take the relative deviation in the corresponding -> 37r 
subtraction constant as an upper limit on what we deem 
acceptable for lu — > 3tt. The sum-rule values for b, see 
Eq. (24), in a; ^ Stt are 



= {0.54( 



0.56 e' 



0.13i 



} 



(43) 



for the Bern and Madrid tttt phase solutions, respectively. 
Note that the phase aigb is significantly smaller for the 
td decay compared to the 0: the imaginary part in the 
subtraction constant is a three-particle-cut effect and as 
such proportional to the phase space available for the three 
pions in the corresponding decays. Enlarging b by the same 
factors as required in the fits to the KLOE (j> -^ Sir data, 
compare Table 4, we evaluate the twice-subtracted u — > 
37r amplitude with 



6={0.83eOO9\0.83eOO7^} 



(44) 



instead. The results for the Dalitz plot parameters tend to 
lie only slightly above the ranges of variation found within 
the uncertainty band for the once-subtracted representa- 
tion; they are quoted in brackets for the full solution in 
Table 5. 

We conclude that we need at least two parameters 
to obtain a polynomial parameterization of the w — > 37r 
Dalitz plot at the 1% level accuracy, and thus directional 
information (owing to the (/)-dependence) is required. Let 
us compare this situation to the rj — > 37r^ Dalitz plot, 
which has the same three-fold symmetry and therefore 
a Dalitz plot distribution that is almost flat. The slope 
parameter is now measured to excellent precision to be 
a{T] -^ 37r°) = -0.0315 ± 0.0015 [44]; we predict a{uj ^ 
Stt) to be about 2.5 times as large, and of opposite sign. 
No higher-order Dalitz parameters have ever been deter- 
mined for 7] -^ 37r°; a theoretical calculation predicts 
/3{r] ^ 37r") = (-4.2 ± 0.7) x lO'^, 7(7? -^ 37r°) = (1.3 ± 
0.4) X 10~^ [7], hence these terms beyond the linear term 
in z will modify the Dalitz plot density only at the few- 
permille level (z < 1). In contrast, /3(lo — > Stt) is larger 
than /3(?7 -^ 37r°) by almost an order of magnitude, see 
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Fig. 12. Contribution to the tttt P-wave inelasticity due to ujn 
intermediate states. Each gray blob corresponds to an oj — >■ Svr 
amplitude, analytically continued to the scattering region. 



Table 5, and hence expected to be significantly more im- 
portant/more likely to be determined experimentally. 

As a final remark on potential uncertainties in the 
Dalitz plot parameters, we note that in an experimen- 
tal investigation of a; — >■ Stt, the invariant mass of the 
detected three pions is going to vary within the natural 
width of the w. We have checked that taking this variation 
into account in the calculation of the decay amplitude the 
two leading Dalitz plot parameters a and l3 change at the 
permille level, hence way below the level of uncertainty of 
our prediction, and probably also significantly below the 
accuracy of any experimental determination in the near 
future. Note that M^ is not changed in the definition of 
the Dalitz plot variables, Eq. (40). 



6 TTTT P-wave inelasticity 

It is well known that KK intermediate states play the 
dominant role in the generation of inelastic effects in the 
TTTT isospin / = S partial wave, where the KK thresh- 
old almost coincides with the position of the /o (980) res- 
onance, up to at least 1.3 GeV. In contrast, the contri- 
bution of KK to the TTTT / = 1 P-wave inelasticity rjl 
is almost negligible, and 7]l is believed to be dominated 
by 47r contributions. Phcnomcnologically, one finds that 
the onset of this inelasticity happens roughly at ^/s ~ 
M^ + A/jr ~ 0.922 GeV, such that it is a natural ques- 
tion to ask whether ujtt intermediate states as an effective 
two-body description of four pions give an adequate de- 
scription of the TTTT P-wave inelasticity, at least in a cer- 
tain energy region above threshold. This requires knowl- 
edge of the inelastic scattering amplitude tttt — )• cjtt in the 
P wave — precisely the crossed process of w — >■ Stt. The 
inelasticity parameter tjI can be calculated from the cut 
contribution shown in Fig. 12, yielding 



; 0.90 - 



0.85- 
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Fig. 13. Inelasticity in the tttt P wave as a function of -^s. The 
dark gray band represents the contribution of ujn intermediate 
states (within our theoretical uncertainty); the light gray band 
shows the error range of the phenomenological determination 
of Ref. [21], with the central value and error limits given by 
the dotted curves. The dash-dotted curve shows the inelasticity 
contribution of the KK intermediate state for comparison. The 
vertical dashed line denotes the cutv threshold. 



seems to provide most of the inelasticity given as the cen- 
tral value of the analysis in Ref. [21]. There is seemingly 
a difference in the threshold behavior, with the ujtt contri- 
bution rising more slowly near threshold; this is due to the 
fact that our tttt — > ujtt partial wave /i(s) happens to have 
almost a zero around 1 GeV. We hasten to add that lutt 
is not the only way to cluster four pions into an effective 
two-body P-wave state; other possibilities like pa would 
just be expected to set in at an even higher effective mass. 
For comparison, we also plot the inelasticity contribu- 
tion from KK intermediate states, derived from the pa- 
rameterization of the TTTT — >■ KK partial wave gl{s) given 
in Ref. [50] (see there for definitions), based on the data 
from Ref. [51]. This yields a tttt P-wave inelasticity con- 
tribution according to 

(46) 
which, in Fig. 13, is indeed seen to be very small, and 
actually remains so up to at least ^/s ~ 2 GeV. 



<llAs)<lL 



i^\f,{sf6{s-{M^+A,Uy) 



(45) 
where ql.is) = X{s, Ml M^)/As, and /i(s) = T{.s)+fi.s) 
is the P-wave projection of the tttt — >■ ujtt amplitude. 

The resulting inelasticity is shown in Fig. 13, up to 
Y^ =1.3 GeV, and for the uncertainty band in the w — >■ 37r 
amplitude discussed in Sect. 4.4. It is seen that this error 
grows rapidly with energy. Phenomenological determina- 
tions of rjl{s) are also plagued by rather large uncertain- 
ties [21,22]; Fig. 13 shows that the ujtt intermediate state 



7 Summary and conclusions 

In this article we have performed a dispersive analysis of 
the decays w — >■ 37r and (f) — > 3tt. Our framework allows 
for a treatment of crossed-channel two-body rescattering 
effects fully consistent with analyticity and unitarity. It 
contains the three-particle cuts generated by My > 3Mt, 
thus going fundamentally beyond an isobar-like descrip- 
tion. We have shown that crossed-channel rescattering 
produces a significant effect on the Dalitz plot clearly ex- 
ceeding uncertainties generated by the phenomenological 
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input for the P-wave phase shift, discontinuities in higher 
partial waves, and elastic P-wave resonances other than 
the p{770). 

Comparing to the very precise measurements of the 
— >■ Stt Dalitz plot by the KLOE collaboration, we have 
found indications for these rescattering effects. Indeed, we 
obtain a considerably improved x^ with the full dispersion 
calculation with one subtraction compared to an isobar- 
like description by Omnes functions. Performing a fit of a 
more flexible, yet less predictive twice-subtracted disper- 
sion calculation, we can achieve an excellent description 
with high probability of both the KLOE and the CMD-2 
data sets, with the more precise KLOE data suggesting 
a significant deviation of the second subtraction constant 
from its sum-rule value. 

For w — )■ Stt we give predictions for the Dalitz plot 
parameters for the full dispersion calculation as well as 
neglecting crossed-channel effects. A precise measurement 
of these parameters that is to be expected from upcoming 
experimental analyses by KLOE and WASA should be 
able to distinguish between both approaches. Moreover we 
have observed that a polynomial fit to the cj — > Stt Dalitz 
plot at percent-level precision should require at least two 
parameters, in contrast e.g. to the decay 77 — > Stt^ that 
displays a similarly symmetric Dalitz plot. The sign of 
the leading polynomial term is constrained by the p peak 
to be positive. 

The inelasticity of the tttt P wave is dominated by Air 
intermediate states, which can be effectively approximated 
by a two-body description as tun. We have found that the 
inelasticity obtained from the analytic continuation of our 
uj -^ Sir amplitude is in reasonable agreement with the not 
very well-constrained phenomenological determination. 

Finally, we wish to remark that it has become amply 
clear in the course of this article that the treatment of both 
decays, uj —^ 3tt and ^ Stt, runs strictly in parallel and 
shows no formal difference other than the mass of the de- 
caying vector meson, and the overall normalization of the 
Dalitz plot, which relates to the (partial) decay widths. 
Obviously, the very same formalism can also be applied 
to the Dalitz plot description of e+e^ — > Stt at arbi- 
trary invariant mass of the electron-positron pair ^s^+e- : 
which replaces the vector mass My. It only has to be 
kept in mind that this would involve an Sg+g- -dependent 
normalization — dispersion theory does not allow us to pre- 
dict the energy-dependence of the e+e~ — ?• Stt total cross 
section. Given that we have shown that crossed-channel 
rescattering effects differ in a non-trivial way for the w 
and (p decays, it may still be interesting to investigate 
the Se+e- -dependence of these corrections in a system- 
atic way. Ultimately, for such a task it might be useful to 
strive for a combination of dispersion theory with micro- 
scopic models that can provide parameterizations for the 
energy-dependent subtraction constants. 
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A On the size of higher partial waves 

In this Appendix, we investigate the potential uncertainty 
in our amplitude solution due to the omission of disconti- 
nuities in higher partial waves for a; — > Stt, the next higher 
one being the F wave. We wish to emphasize to begin with 
that higher partial waves do not vanish for the solution 
J^(s,t, m): projecting onto the F wave (in the s-channel) 

h{s) = -^j dzs{bzt-(izl + l)F{s,t,u) , (47) 

where Zs — cos 6s = {t — u)/K{s), yields contributions from 
the t- and u-channel P-wave amplitudes. /3(s) thus calcu- 
lated just happens to be real. To generalize our approach 
and also include F-wave discontinuities, we may amend 
the decomposition (9) according to 

F{s, i, u) = F(s) + F{t) + F{u) 

+ p!i{zs)g{s) + p^{zt)g{t) + P3{zu)g{u) , (48) 

where 5(s) again only has a right-hand cut, and zt ~ 
(s — u)/K(t), Zu = [s — t)/ k{u). The discontinuities of P 
and F wave are expressed by the relations 

disc/i(s) = discJ^(s) 

= 2i{F{s) + P{s))e{s - 4M^) smSlis)e-'^'^'''> , 
disc /a (s) = disctj(s) 

= 2i{g{s) + gis))0{s - 4M2) sin(5i(s)e-'^3(^) , 

(49) 

where (5|(s) is the tttt F-wave phase shift, and the inho- 
mogeneities are now given by 

:Fis) = mi-zf)i^ + p;g))is), 

gis) ^ -l{{5zt - &zl + 1)(J- + P^g))(.) . (50) 

Note that, in Eq. (50), the notation of angular averaging 
(. . .) is generalized compared to Eq. (16) in the sense that 
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the argument of Pg is taken to be zt, reexpressed in terms 
of s and Zg. We will now proceed to estimate the poten- 
tial size of the F-wave discontinuity contribution G{s) by 
calculating the effects of the p3(1690) resonance. 



A.l The p3(1690) 

The spin-3 resonance p3(1690) is described in terms of a 
totally symmetric third-rank tensor field p^^\ — p't^\T'^ 
(of isospin / = 1), which is subject to the constraints 



d^p^^x = , g'^V^.A = 



(51) 



The numerator of its propagator in momentum space in- 
volves the polarization sum [53] 

pl^l 1/21^3 _ \ ^ J Vi^l Yi'2Vl^3 (y Yl'l'^2 Vl^S 

^^^2^3 - 6 Z^ i^A'l^A'2^M3 ^ 5(^-^MlA'2^ -^M3 



^^"1 Y Y'^2'^3 I Y Y^2 yii's 



"^afi 



lap- 



'-/il^'-P2M3" 
^"P3 



(52) 



where the sum runs over all possible permutations of 
the indices {vi, 1^21 1^3} ■ The simplest p^ interaction La- 
grangians with tttt and nuj are given by 



r — 

'-Pa ~ 



<?/93 7r7r 



{p^.x[d^7r,d-'d^ir]) 



^9f^^XaMl^p^^^Q^Q^^)Q^^Q^^^^ (53) 

where tt = tt^t" denotes the isotriplet of pion fields, and 
uj^ the isosinglet w vector field. From Eq. (53), we can 
calculate the partial decay widths of the pa into tttt and 
■Kuj to be 



r{p?, -^ tttt) 
r{,P3 -> ttw) 



9^ 

y pSTTTT 



(n,f2 _AM^V^'^ 



- IT P3 



2\7/2 



134407rF,2^/;^ 



A(^f^3,M^,A42)'^^ (54) 



which show the expected phase-space dependence for the 
decay of a spin-3 particle. From Mp3 = (1688.8±2.1)MeV, 
Fp^ = (161 ± 10)MeV, B[pz -^ tttt) = (23.6 ± 1.3)%, 
B(p3 -^ TTw) = (16 ± 6)%, and F^ = 92.2 MeV [44], we 
therefore deduce 

|5P3--I = 0.056 ± 0.005 , |5p3,^| = (1.2 ± 0.5) GeV-^ . 

(55) 
Having fixed the coupling constants of the Lagran- 
gian (53), we can proceed to calculate the impact of the 
P3 resonance on the decay uj — >■ 37r. The exchange of a p3 
in the s-channel, see Fig. 14, yields a contribution to G{s) 
as defined in Eq. (48) according to 



Qp^ (s) = Cf 



Ml -s M^ ' 



C, 



5p3 7r7r3p3 7r[^-'"(^ 



60F3M23 



(56) 




Fig. 14. Feynman diagram for the p3(1690)-exchange contri- 
bution to UJ — >■ Stt (in the s-channel). 



K^{s)/M^ is a dimensionless kinematic factor characteris- 
tic for the F wave. The remaining effective coupling Cp is 
numerically evaluated from the above to be 



\C. 



(1.4 ±0.5) X IQ-'^F-^ 



(57) 



For illustration, we compare this to the simplified P-wave 
amplitude J-{s) as given by vector-meson dominance, e.g. 
according to the hidden-local-symmetry (HLS) formal- 
ism [54,55], 



•^HLs(s) = Cp 



Ml 
A'Pp - 



Cp^ 



Nc9 



(58) 



where A'^c = 3 is the number of colors, g ~ 5.8 is the 
universal vector coupling, and we have used the simplest 
choice for the anomalous HLS couplings, ci — C2 = C3 = 
1 [17,18,55]. In this case, the effective coupling constant 
is 

Cp « 0.22 p-3 . (59) 

Comparing Eqs. (57) and (59), we see that the coupling 
constant alone, i.e. before application of the F-wave phase- 
space factors, suppresses the F-wave vs. the P-wave con- 
tribution by three orders of magnitude. 



A. 2 The 



Stt F wave 



We therefore have strong indication that the w — ?► 37r 
F wave /^(s) = Gis) + Gis) is dominated by the term in 
G{s) given by the projection of the crossed-channel P-wave 
terms. With the simplified amplitude (58), we can calcu- 
late this contribution even analytically, eschewing all com- 
plications of path deformation in the complete description. 
The result is 

2M2 

GHLsis) = -Cp " 



7 fb-QK^ + K^^ 1 + K 15-13^2 



2k 



■log 



1-K 



3 



k(s) 



2A/2 - 3so + s 



(60) 



In Fig. 15, we plot the function Cp^Af^K^2(s) x ^hls(s), 
which is seen to vary between 0.53 and 0.14 in the physical 
region 2M^ < ^/s < M^ - M^. In view of Eq. (56), this 
is to be compared to the scale Cp /Cp « lO^'^ for the p^- 
induced F-wave contribution. The conclusion is that we 



F. Niecknig, B. Kubis and S. P. Schneider: Dispersive analysis of a; — > Stt and (^ — ^ Sir decays 



19 




/s [GeV] 

Fig. 15. Tlie function Cp M^k^'^{s) x^hls(s), characterizing 
the F-wave projection of the t- and u-channel P-wave ampli- 
tudes; for details, see text. The vertical dashed lines indicate 
the boundaries of the physical region for io — > 3n. 



expect neglected terms to yield only percent-level correc- 
tions even to the lj ^ 3tt F wave. 

It remains to ask to vifhat extent the assumption of the 
F wave being purely real is justified. Ref. [22] provides a 
simple analytic parameterization of the tttt F-wave phase 
shift, according to which we find, at the upper limits of 
the phase space accessible in a; — > Stt and </) — > Stt, 

Sl{s ^ {M^ - M^f) ^ 1 X 10-4 , 
<5i(s==(M^-M,)2)«7x 10-4 , 



(61) 



which corresponds to 0.06° and 0.4°, respectively. The F- 
wave phase is tiny in all the accessible phase space; ne- 
glecting the imaginary part is an approximation way be- 
low the accuracy of our decay amplitude representation. 



A. 3 Comparison to the pion— pion F wave 

Finally, we wish to briefly point out that the hierarchy 
suggested for the w ^ Stt F wave — complete dominance 
of t-channel exchange over the s-channel resonance tail at 
low energies — should not come as a big surprise, looking 
at the comparable situation in pion-pion scattering. The 
P3 contribution to the tttt F wave that follows from the 
Lagrangian (53) is given by 






89607rF4 M^^ 



2\3 

T ) 

s 



(62) 



This yields a contribution to the F-wave scattering length 
according to 



(4) 



ypSTTTT 



"' UOttF^M^^ 



2.6 X 10-^M. 



r-6 



(63) 



to be compared with a phenomenological value of aj ~ 
5.4 X 10~^M~^ [20]: the pa contribution to the scattering 



length is suppressed by more than two orders of magni- 
tude. The leading-order prediction in chiral perturbation 
theory at one loop, in contrast, yields [56] 



11 



940807r3F4M2 



0{M^) « 2.0 X IQ-'^M^ 



(64) 



which is entirely given by (t- and u-channcl) loop effects. 



B Properties of the dispersion integral 



B.l The angular integral 



The dispersion relation Eq. (22) is constructed for the 
corresponding elastic scattering process Vtt — > tttt with 
My < 3M^, exploiting unitarity. We will now discuss the 
analytic continuation to the decay process and its impact 
on the way the angular integration in the calculation of 
J^{s) is performed. Firstly, we note that we require the 
single-variable amplitude J^{s) in the complex plane for 
the iteration procedure. Secondly, we need to continue to 
the physical vector-meson masses. That means that we 
have to perforin the following continuation: 



Mv 



lim lim T{s) 



>M, 



phy! 



(65) 



Such a continuation was discussed in Ref. [57], comparing 
the dispersive representation to perturbation theory. Sum- 
marizing this analysis, we have to give My an infinitesimal 
positive imaginary part {My — )• My + ie). Apart from the 
angular average integration Eq. (16), the involved func- 
tions depend trivially on My and hence the continuation 
of these functions is simple. The continuation of the an- 
gular integral on the other hand is non-trivial. 

The first issue can be deduced directly from the defin- 
ing equation of the angular average integrals Eq. (16). 
If the argument of the integrand has a real part greater 
than 4M^, the integration contour may cross the right- 
hand cut of J^(s), necessitating its deformation. This issue 
arises in the continuation to vector-meson masses greater 
than SMt^ as a consequence of the three-particle cut, and 
is therefore an inherent property of a decay (as opposed 
to scattering) amplitude. For closer inspection we rewrite 
Eq. (16) in the following form: 



k{s) J^ 



(s) \ 



k{s) 



(66) 
with the integration limits according to the convention 

M^ -^ M^ + ie 



s±{s, M^ + ie) = s±{s, M^) + ie 



ds±{s,Ml 
dMl 



(67) 



20 



F. Niecknig, B. Kubis and S. P. Schneider: Dispersive analysis of cj — >■ Stt and < 



Stt decays 



'^ Im s±(s) 




•C A Res±(s) 



Fig. 16. Trajectories of the integration limits in Eq. (66); see 
text for details. 



we have (see e.g. Ref. [58]) 

' 3so - s + \k{s)\ + ie , se[4M^,a], 



2s+{s) 



2s_(s) 



3so — s + i|K(s)| , s £ [a,b] , 

. 3so - s- |k(s)| , s e [6,oo] , 



3so — s — \k{s)\ + ie , s G 

3so — s — |k(s)| — ie , s G 
3so - s -z|k(s)| , 



s e 



/r2 Mjr-Mj 



4Ml 



1 
a, 6] , 



-,a 



. 3so 



|«:(5)|, 



s e [6, oo] , 



where 



{My - M^f , 6 = (My + Af„)2 



(68) 



(69) 



Care has to be taken at the singularities of the Jacobian 
(k(s) = 0), as here the above parameterization breaks 
down. 

We now deduce where an integration path deforma- 
tion is needed to avoid crossing the right-hand cut. The 
trajectories of the integration limits in the complex plane 
are depicted in Fig. 16. In the scattering region s S [6, oo] 
(corresponding to [D,E] in Fig. 16) we have the standard 
partial-wave projection with line integrations. In the un- 
physical region s G [a, 6] ([C,D]) we have to deform the 
integration contour in the vicinity of a (C), as Res+(s), 
Res_(s) > 4M^ and Im s+(s) > > Im s_(s). Since 
we have the trajectories of the integration limits at hand, 
we can use these as parameterizations for the integration 
contour and exploit the fact that s+ and s_ meet at D. 
We parameterize the integration from s^i^s) to s^{^s) as 
follows: starting from the point S-.(^s)^ we follow the s_ 
trajectory up to D, then use the s+ trajectory to s^^is). 
In the interval s € {[Ml - Ml) 12, a] ([B,C]), s_{s) is 
below the cut, while s+{s) is above. In this part of the de- 
cay region we deform the integration contour as follows: 
the contour follows s_ up to B, then bypasses the branch 
point 4M^ and moves to the upper integration limit above 
the cut via s+. In the interval s e [4M^, (M^ - M^)/2] 
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Fig. 17. Real (solid line) and imaginary (dashed line) part 
of the integrand as a function of ^/s. The occurrence of the 
singularities at the pseudo-threshold (dashed vertical line) is 
explained in the text. 



([A, B]) we can integrate via lines, as the upper and lower 
integration limits lie on the upper rim of the cut. 



B.2 Analytic structure of the integrand 



In this section we wish to briefly discuss the analytic struc- 
ture of the integrand of the dispersion integral in Eq. (22), 
in particular what effects are generated by the fact that 
the vector particle can decay into the final-state particles. 
The analytic structure of the integrand is fully determined 
by the properties of the function k[s) and its zeros. The 
latter are the two-particle scattering thresholds AM^ and 
[Mv + M^f, and the pseudo-threshold [My - M^f . 

Due to the particular structure of the angular integral 
in Eq. (66) the zeros of k(s) in principle give rise to singu- 
larities of square-root order. However, the physical thresh- 
olds at s = AMI ^-nd s = (Mv + Mtt)^ do not contribute to 
the angular integral, since the integration path is shrunk 
to a point. This is not the case at the pseudo-threshold. 
Due to the cut generated by the two-pion threshold and 
the resulting path deformation the integral gives rise to 
a non-vanishing contribution and thus to the aforemen- 
tioned singularity. As explained in the previous section 
the path deformation adheres to the ability of the vec- 
tor particle to decay into the final-state particles, that is. 
My > SMt^. In Fig. 17 we display the integrand of the 
dispersion integral. 

Notice that our observation is in complete agreement 
with a study on the properties of the non-trivial relativis- 
tic two- loop graph given in Ref. [6] . There the singularity 
is explained with a deformation of the path of the an- 
gular integration to infinity as s approaches the pseudo- 
threshold. This is analogous to the previous argument. 
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C Contribution from higher resonances 

In order to estimate the potential influence of higher P- 
wave resonances (/9'(1450), p"(1700)) on the t^ — > Stt de- 
cays, we follow an approach recently suggested in con- 
nection with the pion vector form factor [52]. There, the 
following form factor representation has been fitted to the 
high-precision data for r" 



TT TT^i^T- from Belle [59]: 



F Is) = — — 

^^' M^p - s - iMpFpis) 



exp 



7se 



M^^,-s-iMp,rp,{s) 

5 s e^'f'^ 
Ml„-s-iMp„rp„{s) 



exp 



exp 



967r2F^ 

srp'A^{s) 

irUyiiMl,) 

_ srpnAT,{s) 



(70) 



where 

A^{s) 

rp'.,p"{s) 



= log 



Ml 
M^p 



8Mi 



crUs) log 



l + cr-^ js) 
1 - cr^(s) 



MpS 3 



96ttFI 

Mp,,p. 



s - 4Mi 



^p\p" 



4Mi 



3/2 



r, 



p',p" 



(71) 



We have omitted a kaon-loop contribution to rp{s) (which 
also affects the real exponential multiplying the p prop- 
agator in Eq. (70)) that is retained in Ref. [52]: it does 
not have a large effect on the form factor, and we want to 
translate the above representation into a single-channel 
Omnes form with only elastic tttt rescattering included; 
Eq. (70) treats the p' and p" as purely elastic resonances 
anyway. Besides, inelasticities in the tttt P wave in general, 
and in these higher P-wave resonances in particular, are 
expected to be dominated by 47r rather than KK. Further- 
more, our form of Eq. (70) is only valid above threshold, 
s > AMI. 

Given the treatment of p', p" as elastic resonances, we 
implement these into our formalism by just modifying the 
TTTT P-wave phase shift accordingly and using the phase of 
the form factor in Eq. (70) as the input for a modified 
Omnes function. More precisely, we employ a phase iden- 
tical to the one of Ref. [21] below ^/s = 1.1 GeV, guide this 
smoothly to the phase of the form factor, which is used up 
to y/s = 1.8 GeV (roughly the kinematic range accessible 
in T^ -^ TT^ir'^h'T and therefore fitted in Ref. [52]), before 
guiding the phase smoothly to tt. This phase, calculated 
for the central parameter values of the fit to the Belle 
data [59] given by Mp, = 1.461 GeV, Fp, = 0.353 GcV, 



M. 



p" 



1.732 GeV, Fp„ = 0.141 GeV, 7 = 0.088, S = 
-0.024, (/)! = 0.6, <j>2 = 0.8 [60], is shown in Fig. 18. 

We refrain from investigating the variation within the 
error ranges on these parameters determined in Ref. [52]: 
we regard our estimate of higher-resonance contributions 
as indicative rather (and, as such, they will turn out to 
be entirely negligible) than as a precision study; which, as 




Fig. 18. Alternative P-wave phase shift 5i{s) including the 
effects of p'(1450) and p"(1700) resonances; see text for details. 
The horizontal dashed line denotes the assumed asymptotic 
phase of n, the vertical dashed lines show the boundaries of 
the physical region for the two decays at ^/s — 2M,r , M^j — M^ , 

Ma,-- M^. 



we will argue below, it cannot be as a matter of principle. 
We would like to point out that the form factor repre- 
sentation (70) allows to lead the phase smoothly to n: the 
higher resonances only induce "wiggles" in the phase. This 
is in contrast to the (weighted) sum of Breit-Wigner or 
Gounaris-Sakurai [61] functions as employed in the data 
fit in Ref. [59], which makes the phase rise by tt for each 
resonance, leading to an asymptotic phase of rather 37r in 
the case at hand. 

It turns out that the effect of using an alternative 
Omnes function, based on the phase shown in Fig. 18, 
as the starting point for our dispersive representation of 
the V -^ in Dalitz plots is totally negligible: it changes 
the Dalitz plot distributions by less than 1%, and is there- 
fore smaller than the other uncertainty sources identified 
in the main text. Figure 18 makes this rather plausible: 
the deviations due to the two high-mass resonances show 
up only significantly above the energy range accessible in 
the CO and (f) decays; furthermore, in the integral, the two 
"wiggles" of opposite sign tend to cancel out. Indeed, cal- 
culating e.g. the sum-rule value for the pion charge radius 
according to Eq. (19), we find 



j,2 \V 

sum(p' ,p")/7r 



0.420 fm^ , 



(72) 



very close to the value 0.415 fm^ found with our "standard 
phase" . 

We need to recall the rather indicative nature of this 
investigation. We employ information on the p' and p" res- 
onances gained experimentally from the pion form factor 
measurement in t~ — > tt^tt'^z^t- [59]; it might be possi- 
ble that they couple particularly strongly to wtt, and thus 
have a more pronounced impact on lu — )■ Stt. However, 
phenomenological information does not seem to support 
this: if we compare the pion form factor fit of Ref. [59] 



with an experimental analysis of e^ 



[62], both 
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employing weighted sums of Breit-Wigner or Gounaris- 
Sakurai resonance propagators as fit functions, the pion 
form factor suggests a (modulus of the) coupling strength 
of the p' relative to the p of the order of 15%, while the wtt 
production measurement suggests a p'-to-p ratio of the or- 
der of 10% (despite the fact that here, the p is obviously 
sub-threshold, and phase space tends to emphasize the 
higher-resonance contributions in the cross section). 

We therefore conclude that the influence of higher- 
mass P-wave resonance states on the V —>■ 3tt decays seems 
to be very small, and well contained in the uncertainty es- 
timate we have performed in the main text. 



D Inelasticity parameter 

In this appendix, we wish to briefly derive Eq. (32), fol- 
lowing Ref. [47]. We introduce an inelasticity parameter 
77(5) in the parameterization of the tttt P-wave amplitude 
(see Eq. (31)) according to 

smd{s)e^'^^''> -> ^(1 - r/(s)e-2»*(^)) . (73) 

Thus the unitarity relation for J-{s) takes the following 
modified form 

disc-F(s) = (^(s)+ J-(s))0(s-4M2)(l-ry(s)e-2«*(^)) . 

(74) 
For the homogeneous case J^(s) = 0, we take the logarithm 

disc log J"(s) = 2id{s) — log 77(5) , (75) 

which leads to the modified Omncs solution 



I2(s) = r2(s)exp<j— / ds-^^ 



l&Ml 



s'{s' — s — ie) 



as)S{s)f2{s) 



(76) 



with 



rj ^^^(s) above the cut, 
^(s) = <( ri^/'^{s) below the cut, (77) 

1 elsewhere, 



and 



-'W^-{£/-'^}^ (.8, 

16M2 

For the full solution of the unitarity relation (74) we use 
the product ansatz J^(s) — i7(s)'0(s) to obtain 

, , J-(s)(e^*(«)-r/(s)e-'*'(«)) , ^ 

diBcV'(^) = "^ []' ' . (79) 

Rewriting ^'(■s) into a dispersion relation finally leads to 
the full solution quoted in Eq. (32). 
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